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Abstract 

Numerical simulations of high-strain-rate and high-temperature deformation of pure metals and alloys require 
realistic plastic constitutive models. Empirical models include the widely used Johnson-Cook model and the semi- 
empirical Steinberg-Cochran-Guinan-Lund model. Physically based models such as the Zerilli-Armstrong model, 
the Mechanical Threshold Stress model, and the Preston-Tonks-Wallace model are also coming into wide use. In this 
paper, we determine the Mechanical Threshold Stress model parameters for various tempers of AISI 4340 steel using 
experimental data from the open literature. We also compare stress-strain curves and Taylor impact test profiles pre- 
dicted by the Mechanical Threshold Stress model with those from the Johnson-Cook model for 4340 steel. Relevant 
temperature- and pressure-dependent shear modulus models, melting temperature models, a specific heat model, and 
an equation of state for 4340 steel are discussed and their parameters are presented. 

1 Introduction 

The present work was motivated by the need to simulate numerically the deformation and fragmentation of a heated 
AISI 4340 steel cylinder loaded by explosive deflagration. Such simulations require a plastic constitutive model that 
is valid over temperatures ranging from 250 K to 1300 K and over strain rates ranging from quasistatic to the order of 
10^ Is. The Mechanical Threshold Stress (MTS) model (Follansbee and Kocks 1 14], Kocks 1 29]) is a physically-based 
model that can be used for the range of temperatures and strain rates of interest in these simulations. In the absence 
of any MTS models specificafly for 4340 steels, an existing MTS model for HY-100 steel (Goto et al. (l^ll7|]) was 
initially explored as a surrogate for 4340 steel. However, the HY-100 model failed to produce results that were in 
agreement with experimental stress-strain data for 4340 steel. This paper attempts to redress that issue by providing 
the MTS parameters for a number of tempers of 4340 steel (classified by their Rockwell C hardness number). The 
MTS model is compared with the Johnson-Cook (JC) model (Johnson and Cook 12^ E3l ) for 4340 steel and the 
relative advantages and disadvantages of these models are discussed. 

The MTS model requires a temperature and pressure dependent elastic shear modulus. We describe a number of 
shear modulus models and the associated melting temperature models. Conversion of plastic work into heat is achieved 
through a specific heat model that takes the transformation from the bcc (a) phase to the fee (7) phase into account. 
The associated Mie-Griineisen equation of state for the pressure is also discussed. 

The organization of this paper is as follows. For completeness we provide brief descriptions of the models used 
in this paper in Section |2] Parameters for the submodels required by the MTS model (for example, the shear mod- 
ulus model) are determined and validated in Section 13 Details of the procedure used to determine the MTS model 
parameters are given in Section|3 Predictions from the MTS model are compared with those from the Johnson-Cook 
model in Section|5] These comparisons include both stress-strain curves and Taylor impact tests. Conclusions and 
final remarks are presented in Section0 

*E-mail: banerjee@eng.utah.edu. Phone: (801) 585-5239 Fax: (801) 585-0039 
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2 Models 



In this section, we describe the form of the MTS plastic flow stress model and the associated submodels for the specific 
heat, melting temperature, shear modulus, and the equation of state that have have used for the computations in this 
paper The submodels are used during the stress update step in elastic -plastic numerical computations at high strain 
rates and high temperatures. The submodels discussed in this paper are: 

1 . Specific Heat: the Lederman-Salamon-Shacklette model. 

2. Melting Temperature: the Steinberg-Cochran-Guinan (SCG) model and the Burakovsky-Preston-Silbar (BPS) 
model. 

3. Sliear Modulus: the Varshni-Chen-Gray model (referred to as the MTS shear modulus model in this paper), the 
Steinberg-Cochran-Guinan (SCG) model, and the Nadal-LePoac (NP) model. 

4. Equation of State: the Mie-Griineisen model. 

More details about the models may be found in the cited references. 

The following notation has been used uniformly in the equations that follow. 



e = 


Strain rate 


= 


Plastic strain 


M = 


Shear modulus 


P = 


Current mass density 


Po = 


Initial mass density 




Compression 


ay = 


Yield Stress 


b = 


Magnitude of the Burgers vector 


kb = 


Boltzmann constant 


P = 


Pressure (positive in compression) 


Cp = 


Specific heat at constant pressure 


Cv — 


Specific heat at constant volume 


T = 


Temperature 


T — 


Melting temperature 



Other symbols that appear in the text are identified following the relevant equations. 



2.1 Mechanical Threshold Stress Model 

The Mechanical Threshold Stress (MTS) model (FoUansbee and Kocks fl4ll . Goto et al. 1 17]) gives the following form 
for the flow stress 

gy(£p,£,r) = fTa + {Si(Ti + g'efTe) ^^^''^'' (1) 

where Cq is the athermal component of mechanical threshold stress, Oi is the intrinsic component of the flow stress due 
to barriers to thermally activated dislocation motion, Ce is the strain hardening component of the flow stress, {Si, Se) 
are strain-rate and temperature dependent scaling factors, and po is the shear modulus at K and ambient pressure. 
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The scaling factors Si and Se have the modified Arrhenius form 

hT 



S, - 



Se 



1 - 



In 



1/9. 



1/9. 



1/Pi 



1/Pe 



(2) 



(3) 



where (goi, goe) are normalized activation energies, (eoi, eoe) are constant reference strain rates, and (qi,Pi, qe,Pe) are 
constants. The strain hardening component of the mechanical threshold stress (ce) is given by a modified Voce law 



where 



d(7e „, . 



eo[l~F{ae)]+eiF{ae) 

floo + oio In e + 020 + aaoT 



'1 = aoi 



ail In e + 021 + asiT 
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(5) 
(6) 
(7) 



ln(- 



tanh a 

tanh(Q;) 

hT 

goesb^f^ip, T) 



In 



(8) 
(9) 



and Oq is the strain hardening rate due to dislocation accumulation, 9i is a saturation hardening rate (usually zero), 
(aoj, aij, a2j, a3j, a) are constants (j = 0, 1), is the saturation stress at zero strain hardening rate, coes is the 
saturation threshold stress for deformation at K, goes is the associated normalized activation energy, and eoes is the 
reference maximum strain rate. Note that the maximum strain rate for which the model is valid is usually limited to 
approximately 10^/s. 



2.2 Adiabatic Heating and Specific Heat Model 

A part of the plastic work done is converted into heat and used to update the temperature. The increase in temperature 
(AT) due to an increment in plastic strain (Acp) is given by the equation 



(10) 



where x is the Taylor-Quinney coefficient, and Cp is the specific heat. The value of the Taylor-Quinney coefficient is 
taken to be 0.9 in all our simulations (see Ravichandran et al. |38] for more details on how x varies with strain and 
strain rate). 

The relation for the dependence of Cp upon temperature that is used in this paper has the form (Lederman et al. 
10) 



Ai+ Bit + Ci\t\-°' if T<Tc 
A2 + B2t + C2 t-"' if T>Tc 



T 

t= 1 

Te 



(11) 



(12) 



where Tc is the critical temperature at which the phase transformation from the a to the 7 phase takes place, and 
Al, A2, Bi, B2,a,a are constants. 
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2.3 Melting Temperature Models 
2.3.1 Steinberg-Cochran-Guinan Model 

The Steinberg-Cochran-Guinan (SCG) melting temperature model (Steinberg et al. i4oll ') is based on a modified Lin- 



demann law and has the form 



2a ( 1 - - 
77 



^2(ro-a^l/3). (13) 



where Tmo is the melt temperature at = 1, a is the coefficient of a first order volume correction to the Griineisen 
gamma (Fo). 

2.3.2 Burakovsky-Preston-Silbar Model 

The Burakovsky-Preston-Silbar (BPS) model is based on dislocation-mediated phase transitions (Burakovsky et al. 
L^). The BPS model has the form 



Vw^VJ ; c=|i + -^Pl (14) 




Tyn{p) = t;„(o) 

where ^ is the compression, /ip is the shear modulus at room temperature and zero pressure, /Iq — d^i/dp is the 
pressure derivative of the shear modulus at zero pressure, A'o is the bulk modulus at room temperature and zero 
pressure, ~ dK/dp is the pressure derivative of the bulk modulus at zero pressure, k is a constant, vws is the 
Wigner-Seitz volume, z is the crystal coordination number, a is a constant, and pc{Tm) is the critical density of 
dislocations. Note that ( in the BPS model is derived from the Murnaghan equation of state with pressure as an input 
and may be different from i] in numerical computations. 

2.4 Shear Modulus Models 
2.4.1 MTS Shear Modulus Model 

The Varshni-Chen-Gray shear modulus model has been used in conjunction with the MTS plasticity models by Chen 
and Gray |11] and Goto et al. frT). Hence, we refer to this model as the MTS shear modulus model. The MTS shear 
modulus model is of the form (Varshni [42], Chen and Gray 1 11 1) 

exp(To/T) - 1 

where po is the shear modulus at K, and D, Tq are material constants. There is no pressure dependence of the shear 
modulus in the MTS shear modulus model. 



2.4.2 Steinberg-Cochran-Guinan Model 

The Steinberg-Guinan (SCG) shear modulus model (Steinberg et al. |40|, Zocher et al. j^) is pressure dependent and 
has the form 

/^(P,?") = /^o + ^^Y7^ + - 300); ri^p/po (17) 

where, po is the shear modulus at the reference state(T = 300 K, p = 0, 77 = 1). When the temperature is above Tm, the 
shear modulus is instantaneously set to zero in this model. 
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2.4.3 Nadal-Le Poac Model 



A modified version of the SCG model has been developed by Nadal and Le Poac I35r that attempts to capture the 
sudden drop in the shear modulus close to the melting temperature in a smooth manner. The Nadal-LePoac (NP) shear 
modulus model has the form 



1 



J{T) 



Mo 



where 



J{f) 1 + exp 



9/i p 

dp if-l^ I 

1 + 1/C 



c 



1 + C/(1-T) 



for T 



T 

rn 



e [0,1 + C], 



(18) 



(19) 



fiQ is the shear modulus at K and ambient pressure, ^ is a material parameter, m is the atomic mass, and / is the 
Lindemann constant. 



2.5 Mie-Griineisen Equation of State 

The hydrostatic pressure (p) is calculated using a temperature-corrected Mie-Griineisen equation of state of the form 
(Zocher et al. Ej, see also Wilkins H, p. 61) 

P = + ^qE (20) 

where Cq is the bulk speed of sound, Tq is the Griineisen's gamma at the reference state, Sa = dUs/dUp is a linear 
Hugoniot slope coefficient. Us is the shock wave velocity. Up is the particle velocity, and E is the internal energy per 
unit reference specific volume. The internal energy is computed using 

E=L fc^dT^^^^^^-^ (21) 

where Vb = l/po is the reference specific volume at the reference temperature To. 



3 Submodel Parameters and Validation 

The accuracy of the yield stress predicted by the MTS model depends on the accuracy of the shear modulus, melting 
temperature, equation of state, and specific heat models. The following discussion shows why we have chosen to use 
a temperature-dependent specific heat model, the BPS melting temperature model, the NP shear modulus model, and 
the Mie-Griineisen equation of state model. The relevant parameters of these models are determined and the models 
are validated against experimental data. 

3.1 Specific Heat Model for 4340 Steel 

The parameters for the specific heat model (eauation ll II were fit with a least squares technique using experimental 
data for iron (Wallace et al. f43\, Shacklette 139]) and AISI 3040 steel (fllS). The variation of specific heat with 
temperature predicted by the model is compared with experimental data in Figure \l\ The transition from the bcc a 
phase to the fee 7 phase is clearly visible in the figure. The constants used in the calculations are shown in Tabled If 
we use a constant (room temperature) specific heat for 4340 steel, there will be an unrealistic increase in temperature 
close to the phase transition which can cause premature melting in numerical simulations (and the associated numerical 
problems). We have chosen to use the temperature-dependent specific heat model to avoid such issues. 
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Figure 1: Comparison of experimental data and model prediction of specific heat for 4340 steel as a function of 
temperature. 



Table 1 : Constants used in specific heat model for 4340 steel. 





Ai 


Bi 


Ci a A2 


B2 


C2 a 


(K) 


(J/kg-K) 


(J/kg-K) 


(J/kg-K) (J/kg-K) 


(J/kg-K) 


(J/kg-K) 


1040 


190.14 


-273.75 


418.30 0.20 465.21 


267.52 


58.16 0.35 



3.2 Melting Temperature Model for 4340 Steel 

For the sake of simplicity, we do not consider a phase change in the melting temperature model and assume that the 
iron crystals remain bcc at all temperatures and pressures. We also assume that iron has the same melting temperature 
as 4340 steel. In Figure|2lwe have plotted experimental data (Burakovsky et al. 101, Williams et al. |45], Yoo et al. 
^3l) for the melting temperature of iron at various pressures. Melting curves predicted by the SCG model (Equation 
I13> and the BPS model (Equation 1 14> are shown as smooth curves on the figure. The BPS model performs better at 
high pressures, but both models are within experimental error below 100 GPa. We have chosen to use the BPS melting 
temperature model because of its larger range of applicability. 

The parameters used for the models are shown in Table|2] An initial density (po) of 7830 kg/nx^ has been used in 
the model calculations. 



Wallace et al. (I960) 
□ ASM (1978) 
o Shacklette (1974) 
Model 




600 1000 1400 1800 

T(K) 



3.3 Shear Modulus Models for 4340 Steel 

Figures |3ja), (b), and (c) show shear moduli predicted by the MTS shear modulus model, the SCG shear modulus 
model, and the NP shear modulus model, respectively. Three values of compression (rj = 0.9, 1.0, 1.1) are considered 
for each model. The pressure-dependent melting temperature has been determined using the BPS model in each case. 
The initial density is taken to be 7830 kg/m"^. The model predictions are compared with experimental data for AlSl 
1010 steel and SAE 304 stainless steel. As the figure shows, both steels behave quite similarly as far as their shear 
moduli are concerned. We assume that 4340 steel also shows a similar dependence of shear modulus on temperature. 

The MTS model does not incorporate any pressure dependence of the shear modulus. The pressure dependence 
observed in Figure |3l a) is due to the pressure dependence of T^. Both the SCG and NP shear modulus models are 
pressure dependent and provide a good fit to the data. Though the SCG model is computationally more efficient than 
and as accurate as the NP model, we have chosen to the NP shear modulus model for subsequent calculations for 4340 
steel because of its smooth transition to zero shear modulus at melt. 

The parameters used in the shear modulus models are shown in Table |3] The parameters for the MTS model have 
been obtained from a least squares fit to the data at a compression of 1. The values of /iq and dfi/dp for the SCG 
model are from Guinan and Steinberg 1 18|. The derivative with respect to temperature has been chosen so as to fit 
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Figure 2: Comparison of experimental data and model predictions of melting temperature for 4340 steel as a function 
of pressure. 

Table 2: Parameters used in melting temperature models for 4340 steel. The SCG model parameters are from Gust 
fT^. The bulk and shear moduli and their derivatives have been obtained from Guinan and Steinberg fis'l- The 
parameters for the BPS model at zero pressure have been obtained from Burakovsky and Preston ||7J, Burakovsky 
et al. 10], and the lattice constant (a) is from Jansen et al. |24]. 
Steinberg-Cochran-Guinan (SCG) model 

TmojK) Fq a 

1793 1.67 1.67 

Burakovsky-Preston-Silbar (BPS) model 
~go(GPa) Kg A^o(GPa) fip k z pcjTm) a A vws CA^) a (A) 
~ 166 5.29 81.9 1.8 1 8 0.78 2.9 1.30 0^2 2.865 



the data at a compression of 1 . The NP shear model parameters /iq and C have also been chosen to fit the data. A 
value of 0.57 for C is suggested by Nadal and Le Poac llssll . However, that value leads to a higher value of fi at high 
temperatures than suggested by the experimental data. 

3.4 Equation of State for 4340 Steel 

The pressure in the steel is calculated using the Mie-Griineisen equation of state (equation I20> assuming a linear 
Hugoniot relation. The Griineisen gamma (Fq) is assumed to be a constant over the regime of interest. The specific 
heat at constant volume is assumed to be the same as that at constant pressure and is calculated using equation (|1 1> . 

Figure 13 compares model predictions with experimental data for iron (Bancroft et al. (ll> McQueen et al. |34], 
Barker and Hollenbach |41), mild steel (Katz et al. |28]), 300 series stainless steels (McQueen et al. 1 34]), and for AISI 
4340 steel (Gust et al. ]20]). The high pressure experimental data are not along isotherms and show the temperature 
increase due to compression. The equation of state provides a reasonable match to the experimental data at compres- 
sions below 1 .2 which is reasonable for the simulations of interest in this paper. Improved equations of state should 
be used for overdriven shocks. 

In the model calculations, the bulk speed of sound (Cq) is 3935 m/s and the Unear Hugoniot slope coefficient (So) 
is 1.578. Both parameters are for iron and have been obtained from Brown et al. f^. The Griineisen gamma value (Fq 
= 1 .69) has been interpolated from values given by Gust et al. ] 2QJ . An initial temperature (Tq) of 300 K and an initial 
density of 7830 kg/m^ have been used in the model calculations. 



• Burakovsky et al. (2000) 
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(a) MTS Shear Model 



(b) SCG Shear Model 



g 80 



3 
T3 



u 40 

s: 

m 





AISI 1010 






□ SAE 304SS 






NP (11 = 0.9) 






NP (11 = 1.0) 






NP (11 = 1.1) 



























0.2 0.4 0.6 O.f 

T/T 



(c) NP Shear Model 



Figure 3: Comparison of experimental data with model predictions of shear modulus for 4340 steel. The experimental 
data are for AISI 1010 steel and SAE 304 stainless steel (Fukuhara and Sanpei at standard pressure. 
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Table 3: Parameters used in shear modulus models for 4340 steel . 
MTS shear modulus model 
^iQ (GPa) D (GPa) Tp (K) 

85.0 10.0 298 

SCG shear modulus model 

AJo(GPa) d^i/dp d^i/dT (GPa/K) 

81.9 1.8 0.0387 

NP shear modulus model 

^0 (GPa) dfi/dp C C rn (amu) 

900 TS 004 O080 55.947 
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Figure 4: Comparison of experimental data with model predictions of equation of state for 4340 steel. 

4 Determination of MTS Model Parameters 

The yield strength of high-strength low-alloy (HSLA) steels such as 4340 steel can vary dramatically depending on the 
heat treatment that it has undergone. This is due to the presence of bcc ferrite-bainite phases along with the dominant 
bcc martensite phase at room temperature. At higher temperatures (below the a-j transition) the phases partially 
transform into the fee austenite and much of the effect of heat treatment is expected to be lost. Beyond the transition 
temperature, the alloy is mostly the fee 7 phase that is expected to behave differently than the lower temperature 
phases. Hence, purely empirical plasticity models have to be recalibrated for different levels of hardness of 4340 steel 
and for different ranges of temperature. 

In the absence of relevant microstructural models for the various tempers of 4340 steel, we assume that there is 
a direct correlation between the Rockwell C hardness of the alloy steel and the yield stress (see the ASM Handbook 
1I21I1 '). We determine the MTS parameters for four tempers of 4340 steel. Empirical relationships are then determined 
that can be used to calculate the parameters of intermediate tempers of 4340 steel via interpolation. 

The experimental data used to determine the MTS model parameters are from the sources shown in Table |3 All 
the data are for materials that have been oil quenched after austenitization. More details can be found in the cited 
references. The 4340 VAR (vacuum arc remelted) steel has a higher fracture toughness than the standard 4340 steel. 
However, both steels have similar yield behavior (Brown et al. |6|]). 

The experimental data are either in the form of true stress versus true strain or shear stress versus average shear 
strain. These curves were digitized manually with care and corrected for distortion. The error in digitization was 
around 1% on average. The shear stress-strain curves were converted into an effective tensile stress-strain curves 
assuming von Mises plasticity (see Goto et al. fldl V The elastic portion of the strain was then subtracted from the total 
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Table 4: Sources of experimental data for 4340 steel. 


Material 


Hardness 


Normalize 


Austenitize 


Tempering 


Reference 






Temp. (C) 


Temp. (C) 


Temp. (C) 




4340 Steel 


Rc 30 








Johnson and Cook [26] 


4340 Steel 


Rc 38 


900 


870 


557 


Larson and Nunes [311 


4340 Steel 


Rc 38 




850 


550 


Lee and Yeh |33J 


4340 VAR Steel 


Rc 45 


900 


845 


425 


Chi et al. [12] 


4340 VAR Steel 


Rc 49 


900 


845 


350 


Chi et al. [12] 



strain to get true stress versus plastic strain curves. The Young's modulus was assumed to be 213 MPa. 

4.1 Determination of cr„ 



'a 



The first step in the determination of the parameters for the MTS models is the estimation of the athermal component of 
the yield stress (ffa)- This parameter is dependent on the Hall-Petch effect and hence on the characteristic martensitic 
packet size. The packet size will vary for various tempers of steel and will depend on the size of the austenite 
crystals after the a-7 phase transition. Since we do not have unambiguous grain sizes and other information needed 
to determine a a, we assume that this constant is independent of temper and has a value of 50 MPa based on the value 
used for HY-100 steel (Goto et al. We have observed that a value of 150 MPa leads to a better fit to the modified 
Arrhenius equation for ai and g^i for the Rc 30 temper However, this value is quite high and probably unphysical 
because of the relatively large grain size at this temper. 

4.2 Determination of (Tj and g^i 

From equation Q, it can be seen that ai can be found if ay and ctq are known and ac is zero. Assuming that (Je is zero 
when the plastic strain is zero, and using equation (|2j, we get the relation 





e 



(22) 



Modified Arrhenius (Fisher) plots based on equation (I22t are used to determine the normalized activation energy (^oO 
and the intrinsic thermally activated portion of the yield stress {(Ji). The parameters pi and qi for iron and steels (based 
on the effect of carbon solute atoms on thermally activated dislocation motion) have been suggested to be 0.5 and 1.5, 
respectively (Kocks et al. 1 30], Goto et al. fl^ l. Alternative values can be obtained depending on the assumed shape of 
the activation energy profile or the obstacle force-distance profile (Cottrell and Bilby [ 13], Caillard and Martin ] 10]). 

We have observed that the values suggested for HY-100 give us a value of the normalized activation energy goj for 
Rc = 30 that is around 40, which is not physical. Instead, we have assumed a rectangular force-distance profile which 
gives us values of pi = 2/3 and qi = \ and reasonable values of goi- We have assumed that the reference strain rate is 
eoz = 10«/s. 

The Fisher plots of the raw data (based on Equation \22V ) are shown as squares in Figures (Sja), (b), (c), and (d). 
Straight line least squares fits to the data are also shown in the figures. For these plots, the shear modulus (/i) has 
been calculated using the NP shear modulus model discussed in Sections [2.4.3l and l3.3l The yield stress at zero plastic 
strain {ay) is the intersection of the stress-plastic strain curve with the stress axis. The value of the Boltzmann constant 
(ki,) is 1.3806503e-23 J/K and the magnitude of the Burgers' vector {h) is assumed to be 2.48e-10 m. The density 
of the material is assumed to be constant with a value of 7830 kg/m'^. The raw data used in these plots are given in 
AppendixIXI 

The spread in the data for Rc 30 (Figure |5ja)) is quite large and a very low R^ value is obtained for the fit. This 
error is partially due to the inclusion of both tension and shear test data (in the form of effective tensile stress) in 
the plot. Note that significantly different yield stresses can be obtained from tension and shear tests (especially at 
large strains) (Johnson and Cook lE^l . Goto et al. (3). However, this difference is small at low strains and is not 
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;ure 5: Fisher plots for the intrinsic component of the MTS model for various tempers of 4340 steel. 



11 



Table 5: Values of Ui and g^i for four tempers of 4340 steel. 



Hardness (i?c) 


<y^ (MPa) 


.90i 


30 


867.6 


3.31 


38 


1474.1 


0.44 


45 


1636.6 


1.05 


49 


1752 


1.26 



expected to affect the intrinsic part of the yield stress much. A more probable cause of the spread is that the range of 
temperatures and strain rates is quite limited. More data at higher strain rates and temperatures are needed to get an 
improved correlation for the Rc 30 temper of 4340 steel. 

Figure|5jb) shows the fit to the Fisher plot data for 4340 steel of hardness Rc 38. The low strain rate data from 
Larson and Nunes |31] are the outliers near the top of the plot. The hardness of this steel was estimated from tables 
given in 12111 based on the heat treatment and could be higher than Rc 38. However, the Larson and Nunes |31] data 
are close to the data from Lee and Yeh 1 33] as can be seen from the plot. A close examination of the high temperature 
data shows that there is a slight effect due to the a to 7 phase transformation at high temperatures. 

The stress-strain data for 4340 steel Rc 45 shows anomalous temperature dependent behavior under quasistatic 
conditions. For instance, the yield stress at 373 K is higher than that at 298 K. The fit to the Fisher plot data for this 
temper of steel is shown in Figure|5lc). The fit to the data can be improved if the value of CTq is assumed to be 150 
MPa and qi is assumed to be equal to 2. However, larger values of ctq can lead to large negative values of Ce at small 
strains - which is unphysical. 

The fit to the Fisher data for the Rc 49 temper is shown in Figure|5jd). The fit is reasonably good. More data at 
high strain rates and high temperatures are needed for both the Rc 45 and the Rc 49 tempers of 4340 steel. 

The values of ai and go; for the four tempers of 4340 are shown in Table|5l The value of 50* for the Rc 38 temper 
is quite low and leads to values of the Arrhenius factor {Si) that are zero for temperatures more than 800 K. In the 
following section, we consider the effect of dividing the Rc 38 data into high and low temperature regions to alleviate 
this problem. 

4.2.1 High temperature values of (T; and g^i 

More data at higher temperatures and high strain rates are required for better characterization of the Rc 30, Rc 45, and 
Rc 49 tempers of 4340 steel. In the absence of high temperature data, we can use data for the Rc 38 temper at high 
temperatures to obtain the estimates of ai and g^i for other tempers. We explore two approaches of determining these 
parameters: 

1. Case 1: Divide the temperature regime into three parts: Tq < 573 K; 573 K < To < 1040 K; Tq > 1040 K. 
The values of ai and go; are calculated for each of these regimes from the Rc 38 data. The values of the two 
parameters for temperatures above 573 K are assumed to be applicable to all the tempers. Note that the choice 
of 573 K for the cut-off temperature is arbitrary and loss of temper is likely to occur at a higher temperature. 

2. Case 2: Divide the temperature regime into two parts: Tq < 1040 K and Tq > 1040 K. In this case, we assume 
that the various tempers retain distinctive properties up to the phase transition temperature. All the tempers are 
assumed to have identical values of ai and g^i above 1040 K. 

Case 1: Three temperature regimes. The low and high temperature Fisher plots for Rc 38 4340 steel are shown in 
Figures|6l(a) and (b), respectively. A comparison of Figures|5jb) and|6ja) shows that the low temperature Fisher plot 
has a distinctly lower slope that the plot that contains all the Rc 38 data. The values of ai and goi at low temperatures 
for the Rc 38 temper are 1266 MPa and 1 .85, respectively. The high temperature plot (Figure|6jb)) shows that the slope 
of the Fisher plot is quite steep between 573 K and 1040 K and decreases slightly after the a to 7 phase transition. 
The values of ai and goi for temperatures between 573 K and 1040 K are 1577.2 MPa and 0.371, respectively. After 
the phase transition at 1040 K, these quantities take values of 896.1 MPa and 0.576, respectively. 
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Figure 6: Fisher plots for the intrinsic component of the MTS model for Rc 38 4340 steel assuming three temperature 
regimes. 



Plots of ai and goi as functions of the Rockwell hardness number (for temperatures below 573 K) are shown 
in Figures 0a) and (b), respectively. These plots show a smooth increase in the value of cr; and a decrease in the 
normalized activation energy (goi) with increasing hardness. The high temperature values of g^i for the Rc 38 give 
reasonable values of Si (non-zero) at temperatures above 800 K. However, the lower temperature (less that 1040 K) 
values of the two parameters give a poor fit to the experimental stress-strain data. This is probably due to the anomalous 
behavior of 4340 steel at 373 K and low strain rates. 



Case 2: Two temperature regimes. The two-regime fits to the Fisher plot data for Rc 38 are shown in Figure [s] 
The values of ct; and goi for the Rc 38 temper (in the a phase) are 1528 MPa and 0.412, respectively, while those for 
the 7 phase are 896 MPa and 0.576, respectively. The fits show a jump in value at 1040 K that is not ideal for Newton 
iterations in a typical elastic-plastic numerical code. We suggest that the 7 phase values of these parameters be used if 
there is any problem with convergence. 

Plots of ai and goi as functions of the Rockwell hardness number (for temperatures below 1040 K) are shown in 
Figures |9ja) and (b), respectively. Straight line fits to the ai and goi versus Rc data can be used to estimate these 
parameters for intermediate tempers of the a phase of 4340 steel. These fits are shown in Figure|9] 

The value of ai increases with increasing hardness. However, the value of go; does not decrease monotonically 
with hardness. More experimental data are needed to determine if the trend of goi is physical. The value of goi for 
the Rc 38 temper appears to be unusually low. However, these values lead to good fit to experimental data for Rc 38 
temper For that reason, we have used the two temperature regime values of ai and 50/ for all subsequent computations 
that use these parameters. 



4.3 Determination of a^es and g^e 



Once estimates have been obtained for ai and g^i, the value of Siai can be calculated for a particular strain rate and 
temperature. From equation Q, we then get 



1 

'Sc 



Mo 



(fy - Cq) — Siai 



(23) 



Equation (I23> can be used to determine the saturation value {aes) of the structural evolution stress (cTe). Given a value 
of acs, equation (|9} can be used to compute a^cs and the corresponding normalized activation energy (goes) from the 
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relation 

ln(ae.) = ln(aoe.)~-^^ln(T^) . (24) 

gOesb'^fJ- \<^Oes J 

The value of aes can be determined either from a plot of de versus the plastic strain or from a plot of the tangent 
modulus 9{ae) versus ae- The value of Se is required before can be calculated. Following Goto et al. 1 16], we 
assume that eoe, ioes, Pe, 9e, and goe take the values 10^ /s, 10^ /s, 2/3, 1, and 1.6, respectively. These values are 
used to calculate Se at various temperatures and strain rates. The values of <Ji and goi used to compute a^. vary with 
hardness for temperatures below 1040 K, and are constant above that temperature as discussed in the previous section. 
Adiabatic heating is assumed for strain rates greater than 500 /s. 

Representative plots of <Je versus the plastic strain are shown in Figure [TUf a) and the corresponding 6 versus (7^ 
plots are shown in Figure [T0l" b) (for the Rc 38 temper; strain rate of 1500 /s). Similar plots for the Rc 49 temper for a 
strain rate of 0.0001 /s are shown in Figures fTOI' c) and (d). The plotted value of the tangent modulus (9) is the mean of 
the tangent moduli at each value of de (except for the end points where a single value is used). The saturation stress 
((Tes) is the value at which becomes constant or 9 is zero. Note that errors in the fitting of ai and goi can cause the 
computed value of Ce to nonzero at zero plastic strain. 

The raw data used to plot the Fisher plots for aoes and goes are given in Appendix 151 These data are plotted in 
Figures FTTT a'). (b), (c), and (d). The straight line fit to the data for the Rc 30 temper is shown in Figure [TT]" a). The fit 
to the a phase for Rc 38 4340 steel is shown in Figure [TTT b). Similar Fisher plots for Rc 45 and Rc 49 4340 steel are 
shown in Figures [TTT c) and (d). 

The correlation between the modified Arrhenius relation and the data is quite poor. Considering the fact that special 
care has been taken to determine the value of (Tcs, the poor fit appears to suggest that the strain dependent part of the 
mechanical threshold stress does not follow an Arrhenius relation. However, we do not have information on the error 
in the experimental data and therefore cannot be confident about such a conclusion. We continue to assume, following 
fl^ for HY-100, that a modified Arrhenius type of temperature and strain rate dependence is acceptable for the strain 
dependent part of the yield stress of 4340 steel. 

Values of troes and goes computed from the Fisher plots are shown in Table |5] The value of the saturation stress 
decreases with increasing hardness while the normalized activation energy (at K) increases with increasing hardness. 
For intermediate tempers a median value of 0.284 is assumed for goes and the mean value of 705.5 MPa is assumed 
for CToes- Sti-aight fine fits to the data, as shown in Figures^Ja) and (b), could also be used to determine the values of 
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Figure 8: Fisher plots for the intrinsic component of the MTS model for the a and 7 phases of Rc 38 4340 steel 
assuming two temperature regimes. 
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Figure 10: Plots used to determine the saturation value (ergs) of the structure evolution stress (ije). 
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Figure 1 1 : Fisher plots for the structure evolution dependent component of the MTS model for the a phase of various 
tempers of 4340 steel. 
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Table 6: Values of does and goes for four tempers of 4340 steel. 



Hardness (Re) 


Croes (MPa) 


90es 


30 


1316.1 


0.088 


38 


1058.4 


0.232 


45 


173.5 


0.336 


49 


274.9 


1.245 




25 30 35 40 45 50 25 30 35 40 45 50 

Hardness (R ) Hardness (R ) 



(a) (Toes = -63.9Rc + 3293.4 (MPa) (b) go^ = 0.01656i?c - 0.405 

Figure 12: Values of croes and goes obtained from the Fisher plots for various tempers of the a phase of 4340 steel. 
The dashed lines show the median values of the parameters. 

goes and troes for intermediate tempers of 4340 steel. 

Fits to the data for temperatures greater than 1040 K give us values of coes and goes for the 7 phase of 4340 steel. 
The values of these parameters at such high temperatures are goes = 0.294 and aoes = 478.36 MPa. 

4.4 Determination of hardening rate 6 

The modified Voce rule for the hardening rate (6) (equation is purely empirical. To determine the temperature 
and strain rate dependence of 6, we plot the variation of 6 versus the normalized structure evolution stress assuming 
hyperbolic tangent dependence of the rate of hardening on the mechanical threshold stress. We assume that a — 3. 

Figures Hjl a). (b), (c), and (d) show some representative plots of the variation of 6 with F := ta.nh{aae/aes)/ tanh(a). 
As the plots show, the value of 9i (the value of 9 at F ~ 1) can be assumed to be zero for most of the data. 

It is observed from Figure lOf a) that there is a strong strain rate dependence of 9 that appears to override the 
expected decrease with increase in temperature for the Rc 30 temper of 4340 steel. It can also been seen that 6 is 
almost constant at 298 K and 0.002/s strain rate reflecting linear hardening. However, the hyperbolic tangent rule 
appears to be a good approximation at higher temperatures and strain rates. 

The plot for Rc 38 4340 steel (Figure Hjl b)) shows a strong temperature dependence of 6 with the hardening rate 
decreasing with increasing temperature. The same behavior is observed for all high strain rate data. However, for the 
data at a strain rate of 0.0002/s, there is an increase in 6 with increasing temperature. Figures [13T c) and (d) also show 
an increase in 6 with temperature. These reflect an anomaly in the constitutive behavior of 4340 steel for relatively 
low temperatures (below 400 K) (Tanimura and Duffy |41 ]) that cannot be modeled continuously using an Arrhenius 
law and needs to be characterized in more detail. 

Fits to the experimental data of the form shown in equation (|6} have been attempted. The resulting values of aoo, 
floi, ao2, and ao3 are plotted as functions of Rc in Figures fCT a). (b), (c), and (d), respectively. The points show the 
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13: Plots used to determine 6*0 as a function of temperature and strain rate, a = 3. 
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values of the constants for individual tempers while the solid line shows the median value. These figures show that the 
constants vary considerably and also change sign for different tempers. On average, the strain rate dependence is small 
for all the tempers but significant. The Rc 30 and Rc 45 data points in Figure ^^d) reflect an increase in hardening 
rate with temperature that is nonphysical at high temperatures. 

Instead of using these fits to the experimental data, we have decided to ignore the strain rate dependence of the 
hardening rate and fit a curve to all the data taking only temperature dependence into account (as shown in Figure FTSt . 
Distinctions have not been made between various tempers of 4340 steel to determine this average hardening rate. 
However, we do divide the data into two parts based on the a-7 phase transition temperature. 

The resulting equations for 6*0 as functions of temperature are 

^ [ 15719 - 10.495 T (MPa) for T < IQAQK 
° ^ 1 7516 - 3.7796 T (MPa) for T > IMQK 

This completes the determination of the parameters for the MTS model. 



5 Comparison of MTS model predictions and experimental data 

The performance of the MTS model for 4340 steel is compared to experimental data in this section. In the figures 
that follow, the MTS predictions are shown as dotted lines while the experimental data are shown as solid lines with 
symbols indicting the conditions of the test. Isothermal conditions have been assumed for strain rates less than 500/s 
and adiabatic heating is assumed to occurs at higher strain rates. 

Figure [TSl a) shows the low strain rate experimental data and the corresponding MTS predictions for the Rc 30 
temper of 4340 steel. Comparisons for moderately high strain rates and high temperatures for the Rc 30 temper are 
shown in Figure [TST b). The model matches the experimental curves quite well for low strain rates (keeping in mind 
the difference between the stress-strain curves in tension and in shear). The high strain rate curves are also accurately 
reproduced though there is some error in the initial hardening modulus for the 650 /s and 735 K case. This error can 
be eliminated if the effect of strain rate is included in the expression for 6*0. The maximum modeling error for this 
temper varies between 5% to 10%. 

We have not used the Rc 32 experimental data to fit the MTS model parameters. As a check of the appropriateness 
of the relation between the parameters and the Rc hardness number, we have plotted the MTS predictions versus the 
experimental data for this temper in Figure[n] Our model predicts a stronger temperature dependence for this temper 
than the experimental data. However, the initial high temperature yield stress is reproduced quite accurately while the 
ultimate tensile stress is reproduced well for the lower temperatures. 

The low strain rate stress-strain curves for Rc 38 4340 steel are shown in Figure^] High strain rate stress-strain 
curves for the Rc 38 temper are shown in Figures [T9l a). (b), and (c). The saturation stress predicted at low strain rates 
is around 20% smaller than the observed values at large strains. The anomaly at 373 K is not modeled accurately by 
the MTS parameters used. On the other hand, the high strain rate data are reproduced quite accurately by the MTS 
model with a modeling error of around 5% for all temperatures. 

Experimental data for the Rc 45 temper are compared with MTS predictions in Figures|20|(a) and (b). The MTS 
model underpredicts the low strain rate yield stress and initial hardening modulus by around 15% for both the 173 
K and 373 K data. The prediction is within 10% for the 298 K data. The anomaly at 373 K is clearly visible for 
the low strain rate plots shown in Figure|20ja). The high strain rate data are reproduced quite accurately for all three 
temperatures and the error is less than 10%. 

Comparisons for the Rc 49 temper are shown in Figures|2](a) and (b). The model predicts the experimental data 
quite accurately for 173 K and 298 K at a strain rate of 0.0001/s. However, the anomalous behavior at 373K is not 
predicted and a modeling error of around 15% is observed for this temperature. For the high strain rate cases shown 
in Figure I2ir b). the initial hardening modulus is under-predicted and saturation is predicted at a lower stress than 
observed. In this case, the modeling error is around 10%. 

The comparisons of the MTS model predictions with experimental data shows that the predictions are all within 
an error of 20% for the range of data examined. If we assume that the standard deviation of the experimental data 
is around 5% (Hanson Jl^l) then the maximum modeling error is around 15% with around a 5% mean. This error 
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Figure 14: Variation of the constants fit to the hardening rate equation for various tempers of 4340 steel. 
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Figure 15; Variation of 9o with temperature. 



is quite acceptable for numerical simulations, provided the simulations are conducted within the range of conditions 
used to fit the data. 



6 MTS model predictions over an extended range of conditions 

In this section, we compare the yield stresses predicted for a Rc 40 temper of 4340 steel by the MTS model with those 
predicted by the Johnson-Cook (JC) model. A large range of strain rates and temperatures is explored. In the plots 
shown below, the yield stress (ay) is the Cauchy stress, the plastic strain (e^) is the true plastic strain, the temperatures 
(T) are the initial temperatures and the strain rates e are the nominal strain rates. The effect of pressure on the density 
and melting temperature has been ignored in the MTS calculations presented in this section. The Johnson-Cook model 
and relevant parameters are discussed in Appendixlcl 

6.1 Yield stress versus plastic strain 

Figures l22l' a') and (b) show the yield stress-plastic strain curves predicted by the MTS and JC models, respectively. The 
initial temperature is 600 K and adiabatic heating is assumed for strain rates above 500 /s. The strain rate dependence 
of the yield stress is less pronounced for the MTS model than for the JC model. The hardening rate is higher at low 
strain rates for the JC model. The expected rapid increase in the yield stress at strain rates above 1000 /s (Nicholas 
is^) is not predicted by either model. This error is probably due to the limited high rate data used to determine the 
MTS model parameters. 

The temperature dependence of the yield stress for a strain rate of 1000 /s is shown in Figures l23T a) and (b). Both 
models predict similar stress-strain responses as a function of temperature. However, the initial yield stress is higher 
for the MTS model and the initial hardening rate is lower that that predicted by the JC model for initial temperatures 
of 300K and 700 K. For the high temperature data, the MTS model predicts lower yield stresses. 
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Figure 16: Comparison of MTS prediction with experimental data from Johnson and Cook 12| 
of 4340 steel. 



for the i?c 30 temper 



4340 steel Rc 32 



2000 



Oh 

V 1000 



3 





-o- 0.002/s 298 K 
-w 0.002/s 422 K 
0.002/s 589 K 
-B- 0.002/s 644 K 







0.1 

True Plastic Strain 



0.2 



Figure 17; Comparison of MTS prediction with experimental data from Brown et al. ||6J for the Rc 32 temper of 4340 
steel. 
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Figure 18: Comparison of MTS prediction with experimental data from Larson and Nunes 13 IB for the Rc 38 temper 
of 4340 steel at 0.0002/s strain rate. 



6.2 Yield stress versus strain rate 

The strain rate dependence of the yield stress (at a temperature of 600 K) predicted by the MTS and JC models is 
shown in Figures |53Ia) and (b), respectively. The JC model shows a higher amount of strain hardening than the MTS 
model. The strain rate hardening of the MTS model appears to be closer to experimental observations (Nicholas |36]) 
than the JC model. 

The temperature and strain rate dependence of the yield stress at a plastic strain of 0.3 is shown in Figures l25l a') and 
(b). Above the phase transition temperature, the MTS model predicts more strain rate hardening than the JC model. 
However, at 700 K, both models predict quite similar yield stresses. At room temperature, the JC model predicts a 
higher rate of strain rate hardening than the MTS model and is qualitatively closer to experimental observations. 

6.3 Yield stress versus temperature 

The temperature dependence of the yield stress for various plastic strains (at a strain rate of 1000 /s) is shown in 
Figures |26j a) and (b). The sharp change in the value of the yield stress at the phase transition temperature may be 
problematic for Newton methods used in the determination of the plastic strain rate. We suggest that at temperatures 
close to the phase transition temperature, the high temperature parameters should be used in numerical computations. 
The figures show that both the models predict similar rates of temperature dependence of the yield stress. 

The temperature dependence of the yield stress for various strain rates (at a plastic strain of 0.3) is shown in 
Figures I27f a) and (b). In this case, the MTS model predicts at smaller strain rate effect at low temperatures than 
the JC model. The strain rate dependence of the yield stress increases with temperature for the MTS model while 
it decreases with temperature for the JC model. The JC model appears to predict a more realistic behavior because 
the thermal activation energy for dislocation motion is quite low at high temperatures. However, the MTS model fits 
high temperature/high strain rate experimental data better than the JC model and we might be observing the correct 
behavior in the MTS model. 



6.4 Taylor impact tests 

For further confirmation of the effectiveness of the MTS model, we have simulated three-dimensional Taylor impact 
tests using the Uintah code (Banerjee |2]). Details of the code, the algorithm used, and the validation process have 
been discussed elsewhere (Banerjee 

It is well known that the final length of a Taylor impact cylinder scales with the initial velocity. Figure l28l shows 
some experimental data on the final length of cylindrical Taylor impact specimens as a function of initial velocity. We 
are interested in temperatures higher than room temperature. For clarity, we have separated the high temperature tests 
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Figure 19: Comparison of MTS prediction with experimental data from Lee and Yeh L33.] for the 38 temper of 
4340 steel at high strain rates. 
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Figure 20: Comparison of MTS prediction with experimental data from Chi et al. 
steel. 
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12] for the 45 temper of 4340 



from the room temperature tests by adding an initial internal energy component to the initial kinetic energy density. 
We have simulated three Taylor tests at three energy levels (marked with crosses on the plot). 
The four cases that we have simulated have the following initial conditions: 

1. Case 1: Rc = 30; Lq = 25.4 mm; Dq = 7.62 mm; Ua = 208 m/s; To = 298 K; Source Johnson and Cook Q. 

2. Case 2: R^ = 40; Lq = 30.0 mm; Dq = 6.0 mm; C/q = 312 m/s; To = 725 K; Source Gust id. 

3. Case 3: Rc = 40; Lo = 30.0 mm; Dq = 6.0 mm; Uq = 160 m/s; To = 1285 K; Source Gust d. 

4. Case 4: Rc = 40; Lq = 30.0 mm; Dq = 6.0 mm; Uq = 612 m/s; Tq = 725 K; 

The MTS model parameters for the Rc 30 temper of 4340 steel have been given earlier The MTS parameters 
for the Rc 40 temper of 4340 steel can be calculated either using the linear fit for various hardness levels (shown in 
Figure|9} or by a linear interpolation between the Rc 38 and the Rc 45 values. MTS model parameters at temperatures 
above 1040 K take the high temperature values discussed earlier The initial yield stress in the Johnson-Cook model 
is obtained from the Rc-ctq relation given in AppendixICl 

The computed final profiles are compared with the experimental data in Figures l29l a). (b), (c), and (d). 

For the room temperature test (Figure |29ja)), the Johnson-Cook model accurately predicts the final length, the 
mushroom diameter, and the overall profile. The MTS model underestimates the mushroom diameter by 0.25 mm. 
This difference is within experimental variation (see House et al. fz^ ). 
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Figure 2 1 : 
steel. 



Comparison of MTS prediction with experimental data from Chi et al. 1 12] for the Rc 49 temper of 4340 



The simulations at 725 K (Figure l29l' b)) overestimate the final length of the specimen. The legend shows two MTS 
predictions for this case - MTS (1) and MTS (2). MTS (1) uses parameters (Ti and goi that have been obtained using 
the fits shown in Figure |5] MTS (2) used parameters obtained by linear interpolation between the i?c 38 and Rc 45 
values. The MTS (2) simulation predicts a final length that is slightly less than that predicted by the MTS (1) and 
Johnson-Cook models. The mushroom diameter is also slightly larger for the MTS (2) simulation. 

The final length of the specimen for Case 2 is not predicted accurately by either model. We have confirmed that this 
error is not due to discretization (note that volumetric locking does not occur with the explicit Material Point Method 
used in the simulations). Plots of energy and momentum have also shown that both quantities are conserved in these 
simulations. The final mushroom diameter is not provided by 1 19]. However, the author mentions that no fracture 
was observed in the specimen - discounting a smaller final length due to fracture. In the absence of more extensive 
high temperature Taylor impact data it is unclear if the error is within experimental variation or due to a fault with the 
models used. 

The third case (Figure |55Jc)) was simulated at an initial temperature of 1285 K (above the a-^f phase transition 
temperature of iron). The MTS and Johnson-Cook models predict almost exactly the same behavior for this case. 
The final length is overestimated by both the models. Notice that the final lengths shown in Figure|28]at or near this 
temperature and for similar initial velocities vary from 0.65 to 0.75 of the initial length. The simulations predict a 
final length that is approximately 0.77 times the initial length - which is to the higher end of the range of expected 
final lengths. The discrepancy may be because the models do not predict sufficient strain hardening at these high 
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T = 600 K, various strain rates (MTS) 



T = 600 K, various strain rates (JC) 




Figure 22: 
600 K. 



(a) MTS Prediction. (b) JC Prediction. 

Comparison of MTS and JC predictions of yield stress versus plastic strain at various strain rates for Tq 
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(a) MTS Prediction. 



(b) JC Prediction. 



Figure 23: Comparison of MTS and JC predictions of yield stress versus plastic strain at various strain rates for e 
1000/s. 
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(a) MTS Prediction. 
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(b) JC Prediction. 



Figure 24: Comparison of MTS and JC predictions of yield stress versus strain rate at various plastic strains for Tq = 
600 K. 
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(a) MTS Prediction. 



(b) JC Prediction. 



Figure 25: Comparison of MTS and JC predictions of yield stress versus strain rate at various temperatures for 
0.3 . 



29 




Figure 26: Comparison of MTS and JC predictions of yield stress versus temperature at various plastic strains for 
1000/s. 




Figure 27: Comparison of MTS and JC predictions of yield stress versus temperature at various strain rates for e 
0.3 . 
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temperatures. 

In all three cases, the predictions from the MTS and the Johnson-Cook models are nearly identical. To determine if 
any significant difference between the predictions of these models can be observed at higher strain rates, we simulated 
the geometry of Case 2 with a initial velocity of 612 m/s. The resulting profiles predicted by the MTS and the Johnson- 
Cook models are shown in Figure l29l d). In this case, the MTS model predicts a slightly wider mushroom than the 
Johnson-Cook model. The final predicted lengths are almost identical. Interestingly, the amount of strain hardening 
predicted by the MTS model is smaller than that predicted by the Johnson-Cook model (as can be observed from the 
secondary bulge in the cylinder above the mushroom). We conclude that the Johnson-Cook and MTS models presented 
in this paper predict almost identical elastic -plastic behavior in the range of conditions explored. Please note that quite 
different sets of data were used to determine the parameters of these models and hence the similarity of the results 
may indicate the underlying accuracy of the parameters. 

7 Remarks and Conclusions 

We have determined parameters for the Mechanical Threshold Stress model and the Johnson-Cook model for vari- 
ous tempers of 4340 steel. The predictions of the MTS model have been compared with experimental stress-strain 
data. Yield stresses predicted by the Johnson-Cook and the MTS model have been compared for a range of strain 
rates and temperatures. Taylor impact tests have been simulated and the predicted profiles have been compared with 
experimental data. 

Some remarks and conclusions regarding this work are given below. 

1 . The MTS and Johnson-Cook models predict similar stress-strain behaviors over a large range of strain rates and 
temperatures. Noting that the parameters for these models have been obtained from different sets of experimental 
data, the similarity of the results, especially in the Taylor test simulations, is remarkable. We suggest that this 
is an indication of the accuracy of the models and the simulations. However, the Taylor impact tests show 
that both models predict lower strains at high temperatures than experiments suggest. We are in the process of 
determining paramters for the Preston- Tonks- Wallace model (Preston et al. I37il ') to check if the issue is model 
dependent. 

2. The MTS model parameters are considerably easier to obtain than the Johnson-Cook parameters. However, 
the MTS simulations of the Taylor impact tests take approximately 1.5 times longer than the Johnson-Cook 
simulations. This is partly because the shear modulus and melting temperature models are not evaluated in 
the Johnson-Cook model simulations. Also, the MTS model involves more floating point operations than the 
Johnson-Cook model. The Johnson-Cook model is numerically more efficient than the MTS model and is 
preferable for large numerical simulations involving 4340 steel. 

3. The Nadal-LePoac shear modulus model and the Burakovsky-Preston-Silbar melting temperature model involve 
less data fitting and are the suggested models for elastic-plastic simulations over a large range of temperatures 
and strain rates. The specific heat model that we have presented leads to better predictions of the rate of temper- 
ature increase close to the a-7 phase transition of iron. The shear modulus and melt temperature models are also 
valid in the range of strain rates of the order of 10^ /s. The Mie-Griineisen equation of state should probably be 
replaced by a higher order equation of state for extremely high rate processes. 

4. The relations between the Rockwell C hardness and the model parameters that have been presented provide 
reasonable estimates of the parameters. However, more data for the Rc 30, 45, and 49 tempers are needed for 
better estimates for intermediate tempers. There is an anomaly in the strain rate and temperature dependence of 
the yield strength for Rc 50 and higher tempers of 4340 steel. We would suggest that the values for R^ 49 steel 
be used for harder tempers. For tempers below R^ 30, the fits discussed earlier provide reasonable estimates of 
the yield stress. 

5. The strain hardening (Voce) rule in the MTS model may be a major weakness of the model and should be 
replaced with a more physically based approach. The experimental data used to determine the strain hardening 
rate parameters appear to deviate significantly from Voce behavior is some cases. 
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6. The determination of the values of goes and ages involves a Fisher type modified Arrhenius plot. We have 
observed that the experimental data for the Rc 45 and Rc 49 tempers do not tend to reflect an Arrhenius re- 
lationship. More experimental data (and information on the variation of the experimental data) are needed to 
confirm this anomaly. 
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A Fisher plot data for cr^ and goi 

Tables0|8j|9] and^Jshow the Fisher plot data used to calculate goi and ct; for the four tempers of 4340 steel. 



B Fisher plot data for does and goes 

The data used to compute the parameters coes and goes are shown in Tables [Till 121 II 31 and 1 141 

C Johnson- Cook model and parameters 

The Johnson-Cook (JC) model (Johnson and Cook lEsjl l is purely empirical and has the form 



ay{ep,e,T) — ctq 



B 

1 + -{epY 



[1 + Cln(e*)] [1 - (T*)"] (26) 



e (T - Tr) 
e* = -, T*^j- (27) 

Co [J-m — J-r) 

where ctq is the yield stress at zero plastic strain, and {B, C, n, m) are material constants, eo is a reference strain rate, 
and Tr is a reference temperature. 

The value of cto for 4340 steel in the Johnson-Cook model varies with the temper of the steel. We have fit the yield 
stress versus Rc hardness curve for 4340 steel from the ASM handbook jjUl to determine the value of (Tq for various 
tempers. The equation for the fit is 

(To = exp(Aii?c + A2) (MPa) (28) 

where Ai = 0.0355 In(MPa), A2 = 5.5312 In(MPa), and Rc is the Rockwell-C hardness of die steel. The value of 
B/ao = 0.6339 is assumed to be a constant for all tempers. The strain hardening exponent (n) is 0.26 and the strain 
rate dependence parameter (C) is 0.014, for all tempers. The reference strain rate eo is 1 /s. For temperatures less than 
298 K, thermal softening is assumed to be linear and the parameter m takes a value of 1 . Above 298 K and lower than 
1040 K, TO is assumed to be 1 .03, and beyond 1040 K, to is taken as 0.5 (Lee and Yeh |33|)- Ths reference temperature 
(Tr) is 298 K and the melt temperature (Tm) is kept fixed at 1793 K. These parameters provide a reasonable fit to the 
experimental data presented earlier in the context of the MTS model. 
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Figure 28: The ratio of the final length to the initial length of Taylor impact specimens as a function of initial energy 
density. The experimental data are from Gust fl^, Johnson and Cook |25], Jones and Gillis i27ll . and House et al. 
13 • The tests that we have simulated are marked with crosses. 



Table 7: Fisher plot data used to calculate goi and ai for 4340 steel of hardness Rc 30. 



[fcfcr//ib3ln(£o,/e)]i/''' [i^y-<Ja)/^i]P^ T(K) e (/s) ay (MPa) /i (GP"i) 



0.08333 


0.044658 


298 


0.002 


802.577 


79.745 


0.0782424 


0.0432816 


298 


0.009 


768.052 


79.745 


0.0700974 


0.0425718 


298 


0.1 


750.461 


79.745 


0.0619864 


0.0469729 


298 


1.1 


861.843 


79.745 


0.0408444 


0.0447093 


298 


570 


803.874 


79.745 


0.0747152 


0.0435192 


500 


604 


710.861 


72.793 


0.122804 


0.044074 


735 


650 


648.71 


64.706 
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Jure 29: Comparison of MTS and JC predictions of final Taylor specimen profiles with experimental results. 
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Table 8: Fisher plot data used to calculate goi and cTj for 4340 steel of hardness i?c 38. 







T fKI 

± yrs.) 








077 S4Q? 




?58 


000? 


1 134.12 


81.121 


n OQl 1 1 


L/.L/Jt-x-J It- 






10S7 fn 


7Q 74.S 
1 y , 1 


w.wt / yj 


. J U J T- -7 -7 


9QR 


soo 


1 199 ^8 


7Q 745 

1 y , 1 




V.VJ / UU i / 


9Q8 


1 soo 


1 1 S4. Ifi 

1 1 U 


IQ 74S 
1 y . 1 


1 1 78^1^1 




^7^ 


0009 


1 04.8 8fi 


1 1 , L yj^ 


u.uyuu /oo 


u.uDUoy /o 


J / J 




oJ / .Ul / 


/U.Zol 


0.0819713 


0.0531021 


573 


1500 


910.012 


70.281 


0.134713 


0.0420956 


773 


500 


597.559 


63.398 


0.11695 


0.0461642 


773 


2500 


678.83 


63.398 


0.173098 


0.0367619 


973 


1500 


448.348 


56.515 


0.165137 


0.037095 


973 


2500 


453.773 


56.515 


0.237617 


0.0263176 


1173 


1500 


261.902 


49.632 


0.226689 


0.0287519 


1173 


2500 


291.972 


49.632 


0.322911 


0.0199969 


1373 


1500 


170.886 


42.75 


0.30806 


0.0220299 


1373 


2500 


189.782 


42.75 



Table 9: Fisher plot data used to calculate g^i and di for 4340 steel of hardness Rc 45. 



[hT/i^bHn{eoi/e)]yi^ 


[{ay - a,)/i,r 


T(K) 


e(/s) 


ay (MPa) 


H (GPa) 


0.0514817 


0.0645752 


173 


0.0001 


1429.17 


84.046 


0.0214507 


0.0679395 


173 


1000 


1538.34 


84.046 


0.0934632 


0.0611362 


298 


0.0001 


1255.45 


79.745 


0.038943 


0.0683132 


298 


1000 


1473.83 


79.745 


0.120899 


0.062664 


373 


0.0001 


1260.43 


77.164 


0.0503745 


0.0653759 


373 


1000 


1339.85 


77.164 



Table 10: Fisher plot data used to calculate goi and ct^ for 4340 steel of hardness 49. 



[fc(,T/Ai63ln(eo./e)]i/«- 




T(K) 


e(/s) 


ay (MPa) 


/i (GPa) 


0.0514817 


0.0682425 


173 


0.0001 


1548.31 


84.046 


0.0214507 


0.0711498 


173 


1000 


1645.07 


84.046 


0.0934632 


0.0674349 


298 


0.0001 


1446.46 


79.745 


0.038943 


0.0710307 


298 


1000 


1559.63 


79.745 


0.120899 


0.0653628 


373 


0.0001 


1339.46 


77.164 


0.0503745 


0.0694326 


373 


1000 


1461.75 


77.164 
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Table 11: Fisher plot data used to calculate goes and aoes for 4340 steel of hardness Rc 30. 
[kbT/nbHnjioes/e)] ln(ae«) To (K) T, (K) e (/s) tie. (MPa) (GPa) 
iic = 30 



0.075541 


20.986 


298 


298 


0.002 


1300 


79.745 


0.070454 


19.807 


298 


298 


0.009 


400 


79.745 


0.062309 


19.989 


298 


298 


0.1 


480 


79.745 


0.054198 


19.968 


298 


298 


1.1 


470 


79.745 


0.038728 


20.986 


298 


344 


570 


1300 


78.571 


0.064987 


20.125 


500 


532 


604 


550 


71.984 


0.10315 


19.807 


735 


758 


650 


400 


64.126 



Table 12: Fisher plot data used to calculate f/o, ,s and (To, ,, for 4340 steel of hardness R,. 38. 

[A:,r///i/'lii((,i,,,/V)] ln(a,,) T,, (K) T, (K) ( (/s) ff,, (MPa) f, {GPn) 

i?c = 38 



0.07092 


21.129 


258 


258 


0.0002 


1500 


81.121 


0.08333 


20.986 


298 


298 


0.0002 


1300 


79.745 


0.10779 


21.254 


373 


373 


0.0002 


1700 


77.164 


0.036227 


20.212 


298 


320 


500 


600 


79.183 


0.075873 


20.05 


573 


591 


500 


510 


69.826 


0.11153 


19.552 


773 


785 


500 


310 


63.096 


0.037964 


20.618 


298 


371 


1500 


900 


77.886 


0.070669 


20.367 


573 


614 


1500 


700 


69.246 


0.14023 


20.03 


973 


988 


1500 


500 


56.154 


0.098146 


20.125 


773 


815 


2500 


550 


62.342 


0.13342 


20.088 


973 


995 


2500 


530 


55.989 


0.19164 


19.519 


1173 


1185 


1500 


300 


49.282 


0.25895 


18.891 


1373 


1381 


1500 


160 


42.505 


0.18261 


19.163 


1173 


1193 


2500 


210 


49.047 


0.24641 


19.376 


1373 


1388 


2500 


260 


42.289 



Table 13: Fisher plot data used to calculate goes and croes for 4340 steel of hardness Rc 45. 



[kbT/nbHn{eoes/e)] 


In(aes) 


To(K) 




e(/s) 


aes (MPa) 


II (GPa) 


i?c=45 


0.047192 


19.414 


173 


173 


0.0001 


270 


84.046 


0.085675 


17.034 


298 


298 


0.0001 


25 


79.745 


0.11082 


19.715 


373 


373 


0.0001 


365 


77.164 


0.021177 


19.139 


173 


211 


1000 


205 


83.067 


0.034507 


18.683 


298 


327 


1000 


130 


79.004 


0.043234 


18.826 


373 


397 


1000 


150 


76.553 
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Tabic 14: Fisher plot data used to calculate /yo, ,, and (To, ,, lor 4340 steel of hardness 49. 



[hT/fibHn{eoes/e)] 


ln(<Tes) 


Io(K) 


r.(K) 


e(/s) 


aes (MPa) 


At (GPa) 


Rc=49 


0.047192 


19.254 


173 


173 


0.0001 


230 


84.046 


0.085675 


19.376 


298 


298 


0.0001 


260 


79.745 


0.11082 


19.756 


373 


373 


0.0001 


380 


77.164 


0.02075 


18.951 


173 


207 


1000 


170 


83.17 


0.035325 


19.45 


298 


334 


1000 


280 


78.826 


0.043234 


19.45 


373 


397 


1000 


280 


76.553 
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